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We present an event-driven molecular dynamics study of glass formation in two-dimensional binary mixtures 
composed of hard disks and hard ellipses, where both types of particles have the same area. We demonstrate 
that characteristic glass-formation behavior appears upon compression under appropriate conditions in such 
systems. In particular, while a rotational glass transition occurs only for the ellipses, both types of particles 
undergo a kinetic arrest in the translational degrees of freedom at a single density. The translational dynamics 
for the ellipses is found to be faster than that for the disks within the same system, indicating that shape 
anisotropy promotes the translational motion of particles. We further examine the influence of mixture’s 
composition and aspect ratio on the glass formation. For the mixtures with an ellipse aspect ratio of fc = 2, 
both translational and rotational glass transition densities decrease with increasing the disk concentration at 
a similar rate and hence, the two glass transitions remain close to each other at all concentrations investi¬ 
gated. By elevating k, however, the rotational glass transition density diminishes at a faster rate than the 
translational one, leading to the formation of an orientational glass for the ellipses between the two transi¬ 
tions. Our simulations imply that mixtures of particles with different shapes emerge as a promising model 
for probing the role of particle shape in determining the properties of glass-forming liquids. Furthermore, our 
work illustrates the potential of using knowledge concerning the dependence of glass-formation properties on 
mixture’s composition and particle shape to assist in the rational design of amorphous materials. 

PACS numbers: 61.20.Ja, 82.70.Dd, 64.70.P- 


I. INTRODUCTION 

Despite the fact that glasses are ubiquitous in nature 
and in our daily life, the physical mechanism of the glass 
transition remains one of the most fascinating questions 
in condensed matter physicsii^— Glass formation occurs 
in a variety of materials, ranging from van der Waals liq¬ 
uids and polymers to colloids and granular matter. As¬ 
semblies of spheres provide a prototypical model for ex¬ 
ploring the essential aspects of the liquid-glass transition, 
and much of the physics has been derived from the studies 
of spherical particles, e.g., from the investigations for col¬ 
loidal and granular glasses (see a recent review in Ref. 0) . 
Going beyond the spherical particle paradigm, recent 
great interest has been focusing on probing the glass for¬ 
mation in systems composed of anisotropic particlesi^^— 
From the fundamental point of view, shape anisotropy 
is expected to induce novel dynamical phenomena and 
a better understanding of the glass transition can ben¬ 
efit from exploring the translation-orientation interplay 
in such systems. Investigations for glass formation in 
anisotropic particle systems are also of practical impor¬ 
tance since recent theories^^— and experiment s^ ^ have 
demonstrated that particle shape alone offers a promis- 
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ing route for realizing new kinetically arrested states with 
distinctive mechanical properties. 

A system of hard ellipsoidal particles is among the sim¬ 
plest models of anisotropic particles. The phase behavior 
has been explored extensively by various theoretical ap¬ 
proaches and simulations in both two-dimensional hard 
ellipsesi^^— and three-dimensional hard ellipsoids 
Recently, colloidal ellipsoids have been employed to ex¬ 
perimentally explore a number of novel physical phenom¬ 
ena induced by particle shape.— Although the dynam¬ 
ics in monodisperse hard ellipsoidal particles at low and 
moderate densities has been investigated in recent sim- 
ulations^i^i^ much less information is available con¬ 
cerning the glassy dynamics of hard ellipses^^— and 
hard ellipsoids— The ideal glass transitions in hard el¬ 
lipsoids of evolution (i.e., hard uniaxial ellipsoids) have 
theoretically been treated in the framework of molecu¬ 
lar mode-coupling theory (MCT)— Molecular MCT pre¬ 
dicts that an orientational glass, in which kinetic arrest 
appears in the rotational degrees of freedom but the sys¬ 
tem remains ergodic in the translational degrees of free¬ 
dom, can be formed in hard ellipsoids with large elon¬ 
gations. This implies that a system of hard ellipsoidal 
particles with sufficiently large elongations exhibits two 
glass transitions: The first corresponds to the kinetic ar¬ 
rest in the rotational degrees of freedom and the second to 
the translational freezing of particles. It is only very re¬ 
cently that the formation of an orientational glass and the 
two glass transitions have been observed in experiments 
for monolayers of colloidal ellipsoids with large aspect 
ratios— For hard ellipsoids with small aspect ratio. 
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however, a single glass transition is predicted by molec¬ 
ular MCT^ a prediction that is confirmed by recent 
experiment o^^dO and simulationsMoreover, molecu¬ 
lar MCT suggests that the translational glass transition 
point in hard ellipsoids of evolution displays a nonmono¬ 
tonic variation with increasing aspect ratio4^ Such non¬ 
monotonic behavior is supported by the diffusive proper¬ 
ties of monodisperse hard ellipsoids, as shown in recent 
molecular dynamics (MD) simulations4^ Similar non¬ 
monotonic behavior for the jamming transition point has 
been reported in numerical studies of hard ellipsoidal par- 
ticle a^^d^ and soft ellipses In spite of these recent 
advances, understanding of glassy dynamics is far from 
complete even in two-dimensional hard ellipses. 

The present paper focuses on the glass formation in 
binary mixtures composed of hard disks and hard el¬ 
lipses, where both types of particles have the same size, 
as measured by their areas. Such mixtures offer the po¬ 
tential to better understand the influence of composition 
and particle shape on the glass formation, since crystal¬ 
lization is readily avoided by the anisotropic shape of 
the ellipses. The study is also relevant to colloidal ex¬ 
periments since real synthesized colloidal particles have 
polydispersity not only in size but also in shape. Using 
event-driven molecular dynamics (EDMD) simulations, 
we demonstrate that characteristic glass-formation be¬ 
havior appears upon compression in the translational de¬ 
grees of freedom for both types of particles and in the 
rotational degrees of freedom for the ellipses at appro¬ 
priate compositions and aspect ratios. We analyze the 
glassy dynamics for both species within the same system 
and examine the influence of composition and aspect ra¬ 
tio on the glass formation. Both types of particles are 
found to undergo a kinetic arrest in the translational de¬ 
grees of freedom at a single density, as determined from 
the MCT power-law fits, although the translational dy¬ 
namics for the ellipses is faster than that for the disks 
when comparing their dynamics within the same system. 
The rotational glass transition for the ellipses sets in at a 
density close to the transitional one when the aspect ratio 
k of the ellipses is equal to two, a result that is in agree¬ 
ment with recent experiment o^^d° for monolayers of col¬ 
loidal ellipsoids and simulations for hard ellipses We 
further show that altering mixture’s composition signifi¬ 
cantly affects both glass transition densities when the el¬ 
lipse aspect ratio is fixed. Nevertheless, the density range 
between the two glass transitions appears to be indepen¬ 
dent of composition, because both glass transition den¬ 
sities decrease at a similar rate as the disk concentration 
increases. By elevating the aspect ratio k of the ellipses, 
however, the rotational glass transition density dimin¬ 
ishes at a faster rate than the translational one, leading to 
the formation of an orientational glass for the ellipses be¬ 
tween the two transitions, in line with recent experiments 
for monolayers of colloidal ellipsoid o^^'^^ and MCT pre¬ 
dictions!^ Therefore, our simulations imply that mixing 
particles with different shapes provides a promising way 
for investigating the role of particle shape in determin- 
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FIG. 1. Typical snapshot for a mixture of hard disks and 
hard ellipses with the disk concentration x = 0.5 and the 
ellipse aspect ratio fc = 2 at a high area fraction of 0 = 
0.82. Particles are colored according to their orientation, and 
6 denotes the angle between the particle semi-major axis and 
the a;-axis. 


ing the properties of glass-forming liquids. Furthermore, 
our work illustrates the potential of using knowledge con¬ 
cerning the dependence of glass-formation properties on 
mixture’s composition and particle shape to assist in the 
rational design of amorphous materials. 

The remainder of the paper is organized as follows. 
Section II provides details of our model system and sim¬ 
ulations. Section III begins with a demonstration of ap¬ 
pearance of characteristic glassy dynamics upon increas¬ 
ing density in a mixture of hard disks and hard ellipses, 
and then examines in detail the influence of mixture’s 
composition and aspect ratio on the glass formation. Our 
findings are finally summarized in Section IV. 


II. MODEL AND SIMULATION DETAILS 

The model system considered in the present paper is a 
two-dimensional binary mixture, composed of hard disks 
(denoted by D) and hard ellipses (denoted by E). All 
disks have an identical diameter a. The aspect ratio of 
the ellipses is defined as k = a/6, where a and 6 de¬ 
note the semi-major and semi-minor axes, respectively. 
The particle areas for the two components are chosen 
to be equal (i.e., = 4a6) in the present work. Al¬ 

though other choices are also possible to avoid the crys¬ 
tallization, the present one allows us to better under¬ 
stand the influence of particle shape on the glassy dy¬ 
namics, in particular within the same system, since the 
dynamics for different types of particles is compared at 
the same particle area. Other parameters that specify 
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the system include the number concentration of the disks 
[a; = N£,/{Ni:)+Ne), where {a = D or E) is the num¬ 
ber of a particles], and the area fraction (or density) of 
the system [cj) = 7r(A^D(T^/4 -|- NEah)/L^ with L the box 
dimension]. Figure 1 provides an illustrative snapshot for 
the model system with a; = 0.5 and /c = 2 at a high area 
fraction oi <j) = 0.82. 

Although the disks possess an isotropically spherical 
shape, it is more convenient to assign angular velocities 
to all particles, since the system can then be regarded as 
a binary mixture of hard ellipses with k = 1 for one com¬ 
ponent and fc > 1 for another. Hence, both types of parti¬ 
cles move in both the translational and the rotational de¬ 
grees of freedom in our simulations. EDMD simulations 
are performed in a square box under periodic boundary 
conditions. The reader is referred to Refs. [milliifor 
details about implementing EDMD simulations for hard 
ellipses. All particles have the same mass m and the same 
moment of inertia I. We set m = I = 1 for convenience 
since the present paper only focuses on the general trends 
of static and dynamic quantities in variation with (j), x 
or k. The temperature T is irrelevant for a hard-particle 
system and remains constant due to the conservation law 
of the total kinetic energy. We set the temperature as 
ksT = 1 for simplicity. Length and time are reported 
in units of a and y/ma^/ksT. The results reported in 
the present paper use the system with a total particle 
number of A^ = Ne + Ne = 300. We also performed 
simulations for systems with N = 100 and N = 500 for 
selected state points. Finite size effects are found to be 
small for the properties considered using the present sys¬ 
tem size, and the conclusions given in the present paper 
are unaffected by the system size. A starting configura¬ 
tion with desired (j>, x and k is created using an extended 
Lubachevsky-Stillinger compression algorithm<^“— This 
starting configuration is then equilibrated for at least sev¬ 
eral relaxation times (defined in Subsection III A) before 
collecting data. For each state point, two or four inde¬ 
pendent samples are used to improve the statistics. 


III. RESULTS AND DISCUSSION 

This section first illustrates the appearance of typical 
glassy phenomena upon increasing density in mixtures of 
hard disks and hard ellipses. Specially, we analyze the 
structure and the dynamics for a system where both the 
disk concentration and the ellipse aspect ratio are fixed. 
The influence of mixture’s composition and ellipse aspect 
ratio on the glass formation is then examined in detail. 


A. Glassy behavior in a mixture of hard disks and hard 
ellipses 

A hallmark of the glass transition is that the dynamics 
drastically slows down without any evident change in the 
static pair structure as the liquid approaches its glass 
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FIG. 2. Pair correlation function gdr) for (a) the disks and 
(b) the ellipses, and (c) angular correlation function g 2 ,E{r) 
for the ellipses for mixtures of hard disks and hard ellipses 
with X = 0.5 and fc = 2 at various (f) [indicated in (a)]. The 
power law g 2 {r) ~ is included in (c) to indicate the 

absence of any long-range orientational order in the system 
(see text). 

transition. We illustrate the emergence of the typical 
glassy phenomena in mixtures of hard disks and hard 
ellipses by analyzing the density dependence of structure 
and dynamics for a system with the disk concentration 
X = 0.5 and the ellipse aspect ratio k = 2. 

In order to characterize the positional and the orien¬ 
tational order for each species, we analyze the pair cor¬ 
relation function and the angular correlation function, 
defined by 

’“M = 2^rArNjN^ - 1) “ '’'“J " ‘■“■‘D) ' 

( 1 ) 

and 

92 ,ad) = (cos{2[6»„(r) - 6»a(0)]}), (2) 

where Vaj and daj designate the position and orienta¬ 
tion of particle j belonging to species a, and the en¬ 
semble average < • • • > in Eq. (2) is also performed 
over all particle pairs with the distance r. As shown in 
Figs. 2(a) and 2(b), the long-range positional order for 
both species is apparently suppressed even at the highest 
density studied since gad) quickly decays to unity and 
exhibits several peaks only at short distances. It is also 
obvious from the splitting of the second peak in ga (r) at 
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FIG. 3. (a) Self-intermediate scattering function F’s_Q,(g, f) at 
q = 7 for the disks (dashed lines) and the ellipses (solid lines) 
and (b) 3rd order orientational correlation function Z/ 3 ,£:(t) 
for the ellipses for mixtures of hard disks and hard ellipses 
with X = 0.5 and fc = 2 at various (j)- 


high densities that some short-range order develops upon 
increasing density. This feature is more evident for the 
disks since their isotropic shape favors the formation of 
a local hexagonal structure (see Fig. 1). Nevertheless, 
the overall change in the static pair structure is mild as 
compared to the dramatic change in the dynamics. The 
system also lacks any long-range orientational order. The 
disks are oriented randomly at any density due to their 
spherical shape and hence, the angular correlation func¬ 
tion g 2 ,D{r) for the disks is completely featureless (data 
not shown). Figure 2(c) indicates that the orientational 
order of the ellipses increases upon compression. How¬ 
ever, even at the highest density investigated, the angular 
correlation function g 2 ,Eir) for the ellipses decays much 
faster than that in a nematic liquid crystal,— in which the 
angular correlation function is expected to decay slower 
than the power law g 2 {r) ~ Therefore, the mix¬ 

ture remains disordered in both the translational and the 
rotational degrees of freedom. 

In contrast to the weak changes in the static pair struc¬ 
ture, the dynamics of glass-forming liquids varies sub¬ 
stantially on the approach to the glass transition. The 
typical behavior of a glass-forming liquid is convention¬ 
ally described by the dynamic properties such as the re¬ 
laxation time. The translational dynamics is investigated 
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FIG. 4. (a) and (b) as a function of (j> for 

mixtures of hard disks and hard ellipses with x = 0.5 
and k = 2. Solid and dashed lines are MCT hts for 
the ellipses and the disks, respectively. The fitted re¬ 
sults are as follows. In (a): ( 7 , (()c,Te) = (2.95,0.8294) 

and ( 7 , ((ic.Td) = (3.25,0.8301) for q = 4, (7,<)>c,Te) = 
(2.94,0.8281) and (7,0 c,Td) = (3.13,0.8284) for q = 5, 
(7 ,<(>c,Te) = (3.00,0.8285) and {y,(j)c,To) = (3.33,0.8295) 
for q = 6, (7 ,<(>c,Te) = (3.10,0.8285) and (7 ,<?!>c,Td) = 
(3.43,0.8296) for q = 7. In (b): (7 ,<^c.He) = (2.83,0.8281) 
for n = 2, (7, ((ic.He) = (2.98,0.8280) for n = 3, (7 ,<)<c,He) = 
(3.33, 0.8283) for n = 4. 

by the self-intermediate scattering function, 

\ 

^exp{iq- - r„j(0)]} j , 

i=i / 

( 3 ) 

where i = \/~l and q is the wave number, and the rota¬ 
tional dynamics is explored by the nth order orientational 
correlation function, 

L„,a(t) = ^ ^^exp{m[6l„j(t) - 6l„j(0)]}^ , (4) 

where n is a positive integer. Figure 3 displays F's_Q,(q,t) 
at q = 7 (which is close to the first peak of the static 
structure at high densities) for both species and ,E{t) 
at n = 3 for the ellipses at various area fractions </>. 
The translational relaxation dynamics is found to slow 
down greatly with increasing density for both species. 
In particular, F's_Q,(q, t) exhibits a two-step relaxation 
at sufficiently high area fractions. The two-step decay 
is characteristic of glassy dynamics and implies two re- 
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laxation processes, which are commonly denoted as ,5- 
relaxation and a-relaxationi^ The /^-relaxation corre¬ 
sponds to the rattling motion of particles being trapped 
in cages by their nearest neighbors at short times, while 
the a-relaxation reflects the motion of particles escaping 
from the cages at long times. The rotational dynamics of 
the ellipses is likewise typical of a glass-forming liquid, as 
shown in Fig. 3(b). Obviously, similar glassy dynamics 
does not appear in the rotational dynamics of the disks 
since their isotropic shape allows them to rotate without 
hindrance. This implies that the disks always remain er- 
godic in the rotational degrees of freedom. Our results 
thus indicate a coupling of the translational dynamics 
for both species in mixtures of hard disks and hard el¬ 
lipses, although the structure and the dynamics for dif¬ 
ferent types of particles differ drastically in the rotational 
degrees of freedom. 

It is instructive to compare the translational dynamics 
for both species within the same system because such a 
comparison permits a better understanding for the role 
of particle shape in the translational motion of particles. 
Since the particle areas are chosen to be equal for both 
species, the difference in the translational dynamics is 
solely induced by the particle shape. Figure 3(a) indi¬ 
cates that the dynamics of the ellipses displays a mild 
but detectable speed up as compared to that for the 
disks at low densities [e.g., compare different types of 
lines at (/) = 0.7 in Fig. 3(a)]. The same trend holds for 
high-density states and becomes fairly evident as the liq¬ 
uid enters into the dynamically sluggish regime [see the 
curves at high densities in Fig. 3(a)]. Thus, the trans¬ 
lational motion of particles is promoted by the particle 
anisotropic shape in the same system, a trend that is 
more pronounced at higher densities. The same conclu¬ 
sion remains valid for systems with other compositions 
and other aspect ratios, as discussed in Subsections III 
B and III C. 

It should be emphasized that the translational dynam¬ 
ics for both species is qualitatively the same at any den¬ 
sity despite the presence of quantitative differences. In 
particular, a single translational glass transition point 
4’c,Ta is found to occur for both types of particles, where 
4>c,Tc is determined from the MCT power-law fits of the 
translational relaxation times tt,q for species a, 

TT,a ^ i4‘c,Ta, - (5) 

where 4'c,Ta 7 are fitting parameters, and Tx^a is de¬ 
fined as Fs^a{q,t = TT,a) = 0.1. Figure 4(a) indicates 
that the ^-dependence of is well described by Eq. 
(5)i^ Independent fits produce a common glass transi¬ 
tion point (l)c,Tc = 0.8290 ± 0.0005 at q = 7 for both 
species. Using data for other values of q yields similar 
glass transition densities, which are given in the caption 
of Fig. 4. We note, however, that the fitted exponent 7 
depends on q (see the caption of Fig. 4). Thus, the wave 
number q has merely a quantitative influence on the re¬ 
laxation times as long as q is chosen to be not far away 
from the first peak of the static structure, and the fitted 


glass transition point is independent of q. Likewise, we 
can define the rotational relaxation time for the el¬ 
lipses as Ln,E{t = tr,e) = 0.1, and the MCT power law 
provides a good description for the 4> dependence of tr^e 
at high densities, 

tr,e ^ {(t>c,RE - ■ (6) 

A common rotational glass transition point 4'c,Re = 
0.8281 ± 0.0001 is obtained for the ellipses using data 
for various n [Fig. 4(b)]. Again, the fitted exponent 7 
depends on n (see the caption of Fig. 4). Notice that the 
rotational glass transition sets in at a density quite close 
to the translational one in Fig. 4, a result that arises 
because the ellipse aspect ratio is fixed at k = 2 and that 
is in agreement with the analysis from the self-diffusion 
of monodisperse hard ellipses^S and recent measurements 
for monolayers of colloidal ellipsoids with k « 2j ^^d° The 
density range between the two glass transitions is found 
to be nearly independent of mixture’s composition but 
substantially grow with elevating the aspect ratio of the 
ellipses, as discussed later in detail. 

Therefore, we demonstrate that a mixture of hard disks 
and hard ellipses exhibits typical glass-formation behav¬ 
iors in the translational degrees of freedom for both types 
of particles and in the rotational degrees of freedom for 
the ellipses. The influence of composition and aspect ra¬ 
tio on the glass formation is explored next in Subsections 
III B and III C. 


B. Influence of composition 

Mixture’s composition represents a basic parameter 
for a binary system and has been demonstrated to pro¬ 
foundly affect the glassy dynamics in other types of mix¬ 
tures, such as binary hard spheres^^— and binary hard 
ellipses4i For instance, multiple glasses occur even in 
simple binary hard-disk mixtures, as indicated by MCT 
calculations^iS and simulationsHere, we demon¬ 
strate that the composition strongly influences the glassy 
dynamics in mixtures of hard disks and hard ellipses. 

We consider several disk concentrations ranging from 
X = 0.1 to 0.5, when the ellipse aspect ratio is fixed at 
k = 2. For this aspect ratio, both the positional and the 
orientational order can be suppressed even when the disk 
concentration is low (e.g., x = 0 . 1 ), as evidenced by the 
results in Fig. 5, where we present the pair correlation 
functions for both species and the angular correlation 
functions for the ellipses for various compositions x at 
a fixed area fraction of (/> = 0.82. Figures 5(a) and 5(b) 
also reveal that elevating the disk concentration at a fixed 
density leads to a slight enhancement in the positional or¬ 
der [as reflected by the increasing peak values in ga{r)] for 
both species, in accord with expectations. Dependence 
of the orientational order for the ellipses on composition 
is shown to be quite weak [Fig. 5(c)], implying that com¬ 
position does not play a dominant role in determining the 
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FIG. 5. Pair correlation function ga{r) for (a) the disks and 
(b) the ellipses, and (c) angular correlation function g 2 ,E(r) 
for the ellipses for mixtures of hard disks and hard ellipses 
with k = 2 and various x [indicated in (a)] at </> = 0.82. The 
power law g 2 {r) is included in (c) to indicate the 

absence of any long-range orientational order in the system. 


FIG. 6. (a) Self-intermediate scattering function Fs^a:{q,t) at 
q = 7 for the disks (dashed lines) and the ellipses (solid lines) 
and (b) 3rd order orientational correlation function Ls^sit) 
for the ellipses for mixtures of hard disks and hard ellipses 
with k = 2 and various x at <j) = 0.82. 


spatial alignment of the ellipses. As shown below, how¬ 
ever, altering mixture’s composition significantly affects 
the dynamics of the ellipses in both the translational and 
the rotational degrees of freedom. 

Figure 6 displays the self-intermediate scattering func¬ 
tions for both species and the orientational correlation 
functions for the ellipses for various compositions x at 
(j) = 0.82. Since q and n have merely quantitative ef¬ 
fects on the corresponding time correlation functions if 
they are chosen properly, we focus on Fs^a{q,t) at q = 7 
and Ln^E{t) at n = 3 here and in the following. Fig¬ 
ure 6(a) clearly indicates that the long-time translational 
dynamics for both species significantly slows down upon 
increasing the disk concentration at a fixed density. For 
example, Fs^aiq,t) only exhibits a weak two-step relax¬ 
ation behavior at x = 0.1, while the two-step decay be¬ 
comes fairly evident at a; = 0.5. Accordingly, the transla¬ 
tional relaxation time grows by two orders of magnitude 
when the disk concentration is elevated from 0.1 to 0.5 
at (j) = 0.82. Again, the self-intermediate scattering func¬ 
tion for the ellipses decays faster than that for the disks 
in the same system for all compositions studied [compare 
different types of lines for each x in Fig. 6(a)]. Turning 
to the rotational dynamics, the rotational motion of the 
ellipses is found to similarly be influenced greatly by the 
disk concentration. Specifically, an increase in x similarly 
leads to a significant slowing down of the long-time ro¬ 


tational dynamics for the ellipses at a fixed density [Fig. 
6 (b)]. Thus, adding hard disks into a system of hard el¬ 
lipses tends to inhibit the motion of the ellipses in both 
the translational and the rotational degrees of freedom. 

In order to better illustrate the influence of composi¬ 
tion on the dynamics, we present the ^-dependence of 
the relaxation times for various compositions x in Fig. 
7. Both the translational and the rotational relaxation 
times grow with increasing x at fixed densities, a ten¬ 
dency that becomes more pronounced at higher densi¬ 
ties. Figure 7 also illustrates that the relaxation times 
grow dramatically on the approach to the glass transi¬ 
tion and that the high-density data can be well fitted 
by the MCT power laws (see the solid lines in Fig. 7). 
Thus, the MCT glass transition densities can be deter¬ 
mined from the power-law fits. Again, independent fit¬ 
tings yield a common translational glass transition point 
for both types of particles for each composition. 

Prior studies^^— indicate that the jamming transition 
density in a system of pure ellipses with fc = 2 is signifi¬ 
cantly higher than that in a system of pure disks. Like¬ 
wise, our recent simulations^ indicate that the transla¬ 
tional glass transition in binary hard ellipses with k = 2 
occurs at a much higher density than that in systems of 
pure disks. Hence, it is natural to anticipate that elevat¬ 
ing the disk concentration leads to a diminished trans¬ 
lational glass transition density in the mixture of hard 
disks and hard ellipses with fc = 2, a trend that is al- 
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FIG. 7. (a) TT,a and (b) tji,e as a function of 4> for mixtures 
of hard disks and hard ellipses with k = 2 and various x. Solid 
lines are MCT fits. The fitted glass transition densities are 
shown in Fig. 8. 
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FIG. 9. Pair correlation function Qair) for (a) the disks and 
(b) the ellipses, and (c) angular correlation function g 2 ,E{r) 
for the ellipses for mixtures of hard disks and hard ellipses 
with X = 0.5 and various k [indicated in (a)] at </> = 0.8. The 
power law g 2 {r) ~ jg included in (c) to indicate the 

absence of any long-range orientational order in the system. 


sition where the rotational glass transition shifts from 
higher to lower densities than those corresponding to the 
translational glass transition. This observation might be 
due to uncertainties in the fitting procedures and/or arise 
because the fitted results from the MCT power-law anal¬ 
yses depend on the density range employed. At any rate, 
however, our analysis indicates that mixture’s composi¬ 
tion does not significantly affect the density range be¬ 
tween the two glass transitions. By contrast, the rota¬ 
tional glass transition for the ellipses sets in at a lower 
density than the translational one when the ellipses are 
sufficiently elongated, as discussed in Subsection III C. 


FIG. 8. Translational glass transition density 0c,for both 
species and rotational glass transition density 4‘c,Re fo'' ^^e 
ellipses as a function of x for mixtures of hard disks and hard 
ellipses with k = 2. The line is a guide to the eye. 

ready obvious in Figs. 6 and 7. Such a trend is more 
clearly illustrated in Fig. 8, where the rotational glass 
transition density for the ellipses is shown to similarly 
decrease with increasing the disk concentration. Inter¬ 
estingly, both glass transition densities decrease with the 
disk concentration at a similar rate and remain close to 
each other for all compositions studied. We note that 
the results in Fig. 8 seem to show a slight but consistent 
dependence of the glass transition densities on compo- 


C. Influence of aspect ratio 

We now explore how the ellipse aspect ratio influences 
the glass formation in mixtures of hard disks and hard 
ellipses. For comparative purposes, we fix the disk con¬ 
centration at a; = 0.5. In principle, the influence of as¬ 
pect ratio becomes more evident as the ellipses become 
more elongated and hence, it is crucial to investigate the 
glassy dynamics for mixtures where the ellipse aspect ra¬ 
tio is large enough. However, we have found that different 
types of particles exhibit a strong tendency to separate 
from each other at high densities for large aspect ratios 
(e.g., fc > 4), a behavior that complicates the analysis of 
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FIG. 10. (a) Self-intermediate scattering function Fs^a{q,t) at 
q = 7 for the disks (dashed lines) and the ellipses (solid lines) 
and (b) 3rd order orientational correlation function Z/ 3 ,£;(t) 
for the ellipses for mixtures of hard disks and hard ellipses 
with X = 0.5 and various fc at a fixed area fraction. 


glass formation at any rate. On the other hand, if the 
aspect ratio is too small (e.g., k < 1.5), the rotational 
dynamics of the ellipses resembles that of a free rotator, 
which is quite different from the typical glassy dynamics 
and which has been also observed in hard ellipsoids with 
very small aspect ratios42^ The above considerations 
restrict our attention to the mixtures with 2 < k < 3.5, 
where disks and ellipses are well mixed even at the high¬ 
est densities investigated and the ellipses exhibit typical 
glassy dynamics upon increasing density. 

We first discuss the influence of aspect ratio on the 
static structure of the system. Figure 9(a) indicates that 
the static pair structure for the disks is only weakly influ¬ 
enced by the aspect ratio of the ellipses. The positional 
order of the ellipses displays a slightly diminished trend 
as the aspect ratio elevates [Fig. 9(b)], in line with the 
previous finding^S that shape anisotropy tends to weaken 
the spatial correlations of particles in the translational 
degrees of freedom. In contrast, alterations in the as¬ 
pect ratio lead to evident changes in the structure of the 
ellipses; increasing the aspect ratio enhances the orienta¬ 
tional order of the ellipses [see Fig. 9(c)], because the par¬ 
ticle alignment is entropically more favored for systems 
composed of more elongated particles. Notice that the 
ellipses do not form a nematic phase even for the largest 
k studied, since the angular correlation function still de¬ 
cays much faster than that for a nematic liquid crystal. 
Thus, the ellipses undergo a rotational glass transition at 



FIG. 11. (a) T-r.o and (b) th^e as a function of (j) for mix¬ 
tures of hard disks and hard ellipses with a; = 0.5 and various 
k. The solid lines are MCT hts. The htted glass transition 
densities are shown in Fig. 12. 



FIG. 12. Translational glass transition density (j>c,Ta for both 
species and rotational glass transition density (j)c,RE for the 
ellipses as a function of k for mixtures of hard disks and hard 
ellipses with x = 0.5. The lines are a guide to the eye. 

sufhciently high densities. 

Turning to the dynamics. Fig. 10 exhibits the time cor¬ 
relation functions for various aspect ratios A: at a fixed 
density. Since the rotational motion of the ellipses is 
expected to become more hindered for larger aspect ra¬ 
tios, increasing k inevitably leads to a slowing down of 
the rotational dynamics for the ellipses at fixed densities 
[Fig. 10(b)]. In fact, the influence of the aspect ratio 
on the rotational dynamics of the ellipses is signihcant 
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even at low densities [see Fig. 11(b)], in accord with 
the previous analysis on the basis of diffusive properties 
for hard ellipses.— When the aspect ratio k is elevated 
from 2 to 3.5, the translational dynamics also slows down 
for both species at a fixed density and the difference in 
the dynamics between the two species grows within the 
same system. Since the high-density dynamics traces the 
glass transition point, the trend in Fig. 10(a) appears 
to be consistent with the previous studies on the jam¬ 
ming transition in two-dimensional ellipsoidal particle 
systemsdemonstrating that the jamming transition 
density first grows quickly with increasing k, achieves a 
maximum at a moderate aspect ratio {k Ri 1-5), and then 
slowly becomes smaller for larger aspect ratios. 

Figure 11 depicts the (/)-dependence of the relaxation 
times for various aspect ratios k. The influence of the 
ellipse aspect ratio on the translational dynamics is more 
pronounced at higher densities, while the particle shape 
plays a dominant role in determining the rotational dy¬ 
namics of the ellipses even at low densities. Both the 
translational and the rotational relaxation times grow 
drastically on approaching the glass transition. Again, 
we determine the glass transition points from the MCT 
power-law fits in order to make a quantitative analysis 
for both glass transitions. The variations of the fitted 
glass transition densities with the ellipse aspect ratio are 
displayed in Fig. 12. In line with the above analysis 
for the dynamics in Fig. 10, we observe that both the 
translational and rotational glass transition densities di¬ 
minish as the ellipse aspect ratio elevates. Figure 12 fur¬ 
ther reveals another result: The rotational glass transi¬ 
tion density decreases with the aspect ratio at a faster 
rate than the translational one. Consequently, a density 
regime, where an orientational glass for the ellipses forms, 
emerges due to elevations in the ellipse aspect ratio. This 
regime extends with fc, a trend that is in good agreement 
with recent experiment a^^’^^ for monolayers of colloidal 
ellipsoids and MCT predictions— Therefore, the ellipse 
aspect ratio determines the density range where an ori¬ 
entational glass forms for the ellipses in mixtures of hard 
disks and hard ellipses. 


IV. SUMMARY 

In summary, we explore the glassy dynamics in binary 
mixtures composed of hard disks and hard ellipses. The 
area of the ellipses is chosen to be identical to that of 
the disks, a choice that facilitates elucidating the role 
of particle shape in determining the dynamics since the 
dynamics for both species can be compared within the 
same system at the same particle size. Glass formation 
is demonstrated to occur in the translational degrees of 
freedom for both types of particles and in the rotational 
degrees of freedom for the ellipses under appropriate con¬ 
ditions. The translational dynamics for the ellipses is 
found to be faster than that for the disks within the same 
system. Both types of particles display a common trans¬ 


lational MCT glass transition point despite the presence 
of quantitative differences in their dynamics. We assess 
the influence of composition and aspect ratio on the glass 
formation. Our results indicate that both glass transi¬ 
tion densities are strongly affected by mixture’s compo¬ 
sition. Nevertheless, the density range between the two 
glass transitions appears to be independent of composi¬ 
tion, because both glass transition densities decrease at 
a similar rate as the disk concentration increases. The 
ellipse aspect ratio is shown to profoundly impact on the 
glass formation in mixtures of hard disks and hard el¬ 
lipses. The rotational glass transition density diminishes 
at a faster rate than the translational one when the aspect 
ratio is elevated. Therefore, the density range, where an 
orientational glass forms for the ellipses, appears for suf¬ 
ficiently large aspect ratios. Our simulations thus imply 
that mixtures of particles with different shapes emerge 
as a promising model for exploring the role of particle 
shape in determining the properties of glass-forming liq¬ 
uids. Our work also illustrates the potential of using 
knowledge concerning the dependence of glass-formation 
properties on mixture’s composition and particle shape 
to assist in the rational design of amorphous materials. 
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